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Abstract 

We have investigated the antiferromagnetic insulating phase of the Mott-Hubbard insulator V2O3 
by resonant x-ray Bragg diffraction at the vanadium K-edge. Combining the information obtained 
from azimuthal angle scans, linear incoming polarization scans and by fitting collected data to the 
scattering amplitude derived from the established chemical I2/a and magnetic space groups we 
provide evidence of the ordering motif of anapolar moments (which results from parity violation 
coupling to an electromagnetic field). Experimental data (azimuthal dependence and polarization 
analysis) collected at space-group forbidden Bragg reflections are successfully accounted within 
our model in terms of vanadium magnetoelectric multipoles. We demonstrate that resonant x-ray 
diffraction intensities in all space-group forbidden Bragg reflections of the kind {hkl) m with odd 
h are produced by an E1-E2 event. The determined tensorial parameters offer a test for ab-initio 
calculations in this material, that can lead to a deeper and more quantitative understanding of the 
physical properties of V2O3. 

PACS numbers: 71.30.+h,75.50.Ee,78.70.Ck,78.70.Dm 
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In the case of ions located in crystal positions that are not a center of inversion symmetry, 
hybridization will occur between the valence orbitals of that ion, thus enabling the possibility 
of observing electronic transitions to the hybridized states via the mixed dipole-quadrupole 
(E1-E2) channel in resonant x-ray scattering , which is sensitive to the ordering of 
parity-breaking tensorial moments. For example, the ordering of magnetoelectric toroidal 
moments (anapoles) might be observed. An anapole moment characterizes a system that 
does not transform into itself under space inversion. These toroidal moments were initially 
considered in the context of multipolar expansions in nuclear physics Q and are related to 
a distribution of magnetic fields which is quite different from those produced by parity-even 
multipoles, such as dipole or quadrupole moments. The magnetic field distribution of an 
anapole looks like the magnetic field created by a current flowing in a toroidal winding, anc 
the field is completely confined inside the winding. Parity-breaking E1-E2 contributions 4j 
to scattering can be expressed in terms of polar and magnetoelectric tensors that contain 
the anapole operator l| • The study of these contributions to resonant x-ray diffraction is of 
fundamental importance in current developments of the electronic structure of materials with 
complex electronic properties, such as magnetoelectricity, piezoelectricity and ferroelectricity 
that are of potential technological interest [5|. The cross correlation between magnetism and 
ferroelectricity in materials with coexistence of spontaneous magnetization and polarization, 
termed multiferroics, has recently become a subject of great scientific impact. 

Vanadium sesquioxide, V2O3, considered as a Mott-Hubbard metal-insulator j(| system, 



has been the o 
tal 



11 



12, 



13 



jject of intense study from both theoretical lOj and experimen 



14j points of view in the last decades. This compound has an interesting 



phase diagram, with an antiferromagnetic insulator phase (AFI) at low temperatures, and 
a paramagnetic metallic one (PM) above the Neel temperature (Tjv ~ 150 K). The metal- 
insulator transition is accompanied by a strongly first-order structural phase transition in 
which the room-temperature corundum structure (R3c) is modified to monoclinic I2/a. In 
1978 Castellani, Natoli and Ranniger [jj proposed a theoretical model to explain the mag- 
netic structure in the AFI phase from the ordering pattern of the occupation of the t2 9 
(degenerate) orbitals. This model was considered valid until 1999, when resonant x-ray 
diffraction (Paolasini et al. II]), and magnetic dichroism measurements (Park et al. [13I ]) 
demonstrated that the spin of the Vanadium atoms was Sy — 1, whereas the model of 
Castellani predicted Sy = 1/2. 
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The magnetic structure is such that the I-centering cell translation is time inverting. 
Because of the magnetic moments being colinear in the a m -c m plane, normal magnetic peaks 
are observed only at reflections (hkl) m with even h even and odd (k+l). Yet, peaks have been 
measured at odd h, even [h + I), at the resonant prepeak of Vanadium K-edge. They were 
initially interpreted as produced by orbital (time-even) ordering 

of the observations carried out by Lovesey et al. [h]] and Tanaka |8|] demonstrated that the 
resonant Bragg diffraction intensities are produced by the ordering of magnetic (time-odd) 
V multipoles. There has been controversy on the dominant resonant event producing the 
_d intend, havtng been proposed a^arity-even E2-E2 event Q, a parity-odd 
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121 ] . However, analyses 



E1-E2 event [8j , or a combination of them [9|, ]15j . 

In this paper we present measurements at space-group forbidden Bragg reflections (hkl) m 
with odd h in the low temperature monoclinic phase of V2O3. Together with energy profiles 
and azimuthal dependence, we measured the polarization dependence at a fixed azimuthal 
point varying the incident linear polarization with a phase-plate setup. In this way, we 
can obtain the Stokes parameters for the secondary beam as a function of both incident 
and scattered linear polarization angle, without moving the crystal. This method can elu- 
cidate the presence of resonances that are very close in energy, playin g on their relative 
phase shifts [16J. Experiments were carried out at the ID20 beamline [17|] of the ESRF. 
A single crystal of 2.8% Cr-doped (Vi- x Gt x ) 2 3 (x = 0.028) was mounted on the four- 
circle diffractometer with vertical scattering geometry. For the varying linear incoming 
polarization, a diamond phase-plate of thickness 300 /im was inserted into the incident 
beam. The phase plate was operated in half-wave plate mode to rotate the incident po- 
larization into an arbitrary plane Ij], 19] described by Stokes parameters P\ = cos(2rj) 
and P2 = sin(2r)), being 77 the angle between the incident beam electric field and the axis 
perpendicular to the scattering plane, i.e., 77 = corresponds to polarization perpendicular 
to the scattering plane (a polarization). The Poincare- Stokes parameters are defined by 
P x = (|^| 2 - \s w \ 2 )/Po,P2 = 2Re{e*e„)/P ,P 3 = 2lm(s* a e n )/P Q , with P = \e a \ 2 + \e n \ 2 the 
total intensity, and where e a and are the components of the beam polarization vector 
perpendicular and parallel to the scattering plane. 

Fig. 1 shows the energy profiles for the (302) m and (102) m reflections. They contain 
a single resonant peak centered around 5.465 keV, as it is expected for reflections (hkl) m 
with odd h, in which the resonant peak from 5.47 to 5.49 keV, ascribed to an El — El 
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event is forbidden [lOj, Energy profiles in Fig. 1 show the presence of a shoulder at 
E=5.4665 keV for both reflections, which opens the possibility for the existence of two 
lorentzians separated by approximately 2 eV contributing to the observed intensity. This 



energy separation can be related to crystalline electric fielc 



energy transfer. The value of 
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In order to elucidate the 



lODq was estimated as 2.1 eV from RIXS measurements 
possibility of interference of different lorentzians we have performed polarization analysis 
measurements of the dependence of the Stokes parameters of the secondary beam with the 
angle of primary linear polarization 77 using the phase plate (Fig. 2). In the case of the 
(102) m reflection we have repeated the polarization analysis measurements at 3 different 
energies of the incoming x-rays (Fig. 2-b). The presence of multiple lorentzians at different 
energies with different tensorial properties would be revealed by the appearance of circular 
polarization in the polarization measurements, which can be estimated from the dependence 
with linear incident polarization of the measured secondary Stokes parameters P[ and P3, 
and by a variation in the shape of the Stokes parameter curves when the energy of the x-rays 
is changed. We try a fit of the polarization data with a real Jones matrix derived from the 
unit-cell structure factors for the different polarizations F n „> , F^-a 1 , F<j-n', F n _ n i (formulae 



for the outgoing Stokes parameters can be seen in ref. |20| ) together with a depolarization 
factor (1 — l/2<isin 2 77) that multiplies the linear components N\ and N2 of the Stokes vector 
of the incoming light, being d an adjustable parameter. This depolarization produced by 
the phase plate depends on the angle 77 indicating the angle of linear polarization of the 
incoming beam and its maximum is expected at rj = 90 degrees. From the fitting to the 
(302) m data we obtain the following parameters, 

d = 0.162 ±0.011 
F^/Fe-j = 0.8 ± 0.3 
F w _ a , / F a _ a , = 0.64 ±0.02 

F a ^/F a _ a , = -3.01 ± 0.05, (1) 

and from the fitting to the (102) m polarization data we obtain, 

d = 0.362 ±0.014 
F v ^lF a -a> = 0.099 ±0.005 
F^/F^, = 1.20 ±0.01 



F^jF^ = 0.770 ±0.005. (2) 

In both cases the magnitudes of the depolarization factor d is similar to that estimated 
from preliminary measurements with the phase plate. We obtain two different values of the 
depolarization d due to the fact that polarization measurements in the (302) m and in the 
(102) m reflections were performed in different experiments. The good agreement obtained 
in the fittings for both reflections with experimental data for a model in which we take into 
account the depolarization introduced by the phase plate, together with the fact that there is 
no significative dependence of the shape of Stokes parameters curves when the energy of the 
x-rays is changed in the measurements of (102) m reflection leads us to conclude that there is 
no evidence of the appearance of circular polarization and that the measured intensities in 



both reflections are produced by single oscillators with the same tensoria 
this supports the validity of fitting the data with a single oscillator model 



character. All of 
9j. In this aspect, 



the model that will be used to describe experimental data differs to the case of K 2 Cr04, 



where there is a strong evidence of the appearance of circular polarization [16|, [20j and the 
intensities were consequently modelled in terms of different lorentzians centered at different 
energies for the different resonant events. 

In Fig. 3 we show the azimuthal dependence of (302) m and (102) m reflections measured 
at T=100 K in the AFI phase. In order to eliminate the effect of the absorption in the 
azimuthal curves, data have been corrected according to the formula 

^corr = 4bs(l + sina / sinai), 

where ao and a.\ are the incident and reflected angles of the beam with the sample. The 
measured azimuthal dependence of (302) m shows a satisfactory agreement with previously 
published data [9|, as it can be seen in Fig. 3-a. 

The azimuthal dependence and polarization data collected for (102) m reflection present a 
good agreement with the expression for the parity-breaking event resonant x-ray scattering 
F(E1 — E2) presented in ref. 0, being possible to fit together polarization scans curves and 
azimuthal data. The azimuthal and polarization data for the (302) m reflection can also be 
fitted to the El — E2 structure factor expressions together with (102) m data in terms of 
the same set of tensorial parameters, which permits us to conclude that intensities in both 
(302) m and (102) m are produced by an El — E2 event. The result of the fitting is shown in 



5 



Fig. 3. The determined parameters from fitting \F(E1 — E2)\ 2 to the data are, 

lm(G\) /lm(Gl) = 0.263 ± 0.013 

(Gl)/Im(G 3 3 ) = 3.48 ±0.06 
lm(Gl)/Im(G 3 3 ) = -8.38 ± 0.02 
Re(Gl)/lm(Gl) = 3.908 ± 0.013 
Re(G 3 2 ) /lm(Gl) = -0.205 ± 0.007 
lm(Gl) /lm(Gl) = -3.410 ± 0.015 

(G 3 }/lm(Gl) = -4.33 ±0.03. 

(3) 

Values for lm(G\) and (Gj) are direct estimates of the orbital anapolar moment (f2), lm(Gf,), 
Re(Gl), Im{Gf} and (Gq) are estimates of the moment ((L <g> (L <g> ft) 2 ) 3 ), where L is the 
operator for orbital angular momentum, and Im(G|) and R,e(Gf) are estimates of the mo- 
ment ((L (gi n) 2 ), where n is the polar unit vector. Operators O and (L <g> (L <g> f2) 2 ) 3 are 
true spherical tensors while (L <g> n) 2 is a pseudo-spherical-tensor. Additional information 
about parity-odd tensors may be found in refs. l|, y, 21, 22. The presence of a parity-even 



contribution FiE2 — E2) to resonant x-ray scattering in this kind of reflections has been 
suggested {9], [lO, Q and within measured reflections its strongest contribution would appear 
in the reflection (302) m as F(E2 — E2) is weighted by a global multiplicative factor smv, 
with v — hx + ky ± Iz, being x, y, z crystallographic parameters and the h, k, and / the 
Miller indices of the reflection [9j. However, our fitting describes data in (302) m with great 
precision, which makes us conclude that parity-even E2 — E2 contribution is absent to a good 
degree of approximation. From the fitted parameters, we can determine the direction of the 
projection in the a m -c m monoclinic plane of the anapolar moments of the Vanadium ions, 
which from the quotient between lm.(G\) and (Gq) is estimated as forming 4 degrees with 
the (101) m reciprocal lattice vector. In Fig. 4, we show the projection in the a m -c m plane of 
the anapolar moments of the eight vanadium ions in the monoclinic unit cell, together with 



the arrangement of magnetic moments 23j established below T/y. 

In conclusion, our results demonstrate that resonant x-ray diffraction at the V K-edge 
for space-group forbidden Bragg reflections of the kind (hkl) m with odd h is produced by 
a parity-breaking E1-E2 event, being the contribution from parity-even transitions absent 
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to a very good degree of approximation. Polarization analysis measurements at the (302) m 
and (102) m reflections probe the fact that the intensities measured are coming with a good 
degree of approximation from single lorentzians. Experimental data (azimuthal variation 
and polarization analysis) collected at different space-group forbidden Bragg reflections are 
successfully accounted within our model in terms of values for the V magnetoelectric multi- 
poles. The derived tensorial parameters include direct estimates of expectation values of the 
vanadium anapolar moment and other magnetoelectric moments [9j . These results solve the 
controversy on the origin of the resonant x-ray diffraction intensities in V2O3 js, 9, 10, 15]. 
The derived tensorial parameters offer a test for ab-initio calculations of the electronic struc- 
ture of vanadium sesquioxide that can lead to a deeper and more quantitative understanding 
of its electronic properties. 
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FIG. 2: Linear incident polarization dependence of the Stokes parameters P[ and P' 2 of the sec- 
ondary beam in the reflections (302) m (ifi = —39.0 degrees) collected at T=100 K with an incident 
energy of the x-rays E = 5.4660 keV and for the reflection (102) m at if) = —31.1 degrees and T=150 
K and at three different energies E = 5.4655, 5.4660 and 5.4665 keV. Dashed line corresponds to 
the fitting of the data to a general Jones matrix produced by a single oscillator including the effect 
of the phase plate depolarization. Continuous line corresponds to the fitting of polarization data 
together with azimuthal scans using the tensorial parameters presented in eq. (3). In the origin of 
the azimuthal angle the (010) m reflection is in the plane of scattering. 

[18] V. Scagnoli et al., submitted to J. Synchr. Rad. (2009). 

[19] C. Giles, C. Vettier, F. de Bergevin, C. Malgrange, G. Grubel, and F. Grossi, Rev. Sci. 

Instrum. 66, 1518 (1995). 
[20] J. Fernandez- Rodriguez, S. W. Lovesey, and J. A. Blanco, Physical Review B 77, 094441 

(2008). 

[21] P. Carra and R. Benoist, Phys. Rev. B 62, R7703 (2000). 

[22] P. Carra, A. Jerez, and I. Marri, Phys. Rev. B 67, 045111 (2003). 

[23] R. M. Moon, Phys. Rev. Lett. 25, 527 (1970). 



9 




FIG. 3: Azimuthal angle scans measured on the reflections (302) m and (102) m at T = 100 K 
with energy E = 5.46627 keV. In the case of (302) m we show the agreement of new data with 
previously published data [9|]. The continuous line correspond to the fitting to the expressions 
for the F(E1 — E2) scattering length presented in [9|] with the tensorial parameters shown in eq. 
(3). At the origin of the azimuthal angle (010) m reflection is in the plane of scattering. Both 
polarization channels a — a' and a — a' are shown. 
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FIG. 4: (color online) Positions of the Vanadium ions in the monoclinic unit cell adopted by V2O3 
below the Neel temperature, together with the configuration of the magnetic moments 23[ (red 
arrows) and the determined projection of their anapolar moments (black arrows) in the a m -c m 
plane. The basis vector h m is normal to the plane of the diagram. 
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